Singularities in fidelity surfaces for quantum phase transitions: a geometric perspective 
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The fidelity per site between two ground states of a quantum lattice system corresponding to different values 
of the control parameter defines a surface embedded in a Euclidean space. The Gaussian curvature naturally 
quantifies quantum fluctuations that destroy orders at transition points. It turns out that quantum fluctuations 
wildly distort the fidelity surface near the transition points, at which the Gaussian curvature is singular in the 
thermodynamic limit. As a concrete example, the one-dimensional quantum Ising model in a transverse field 
is analyzed. We also perform a finite size scaling analysis for the transverse Ising model of finite sizes. The 
scaling behavior for the Gaussian curvature is numerically checked and the correlation length critical exponent 
is extracted, which is consistent with the conformal invariance at the critical point. 

PACS numbers: 03.67.-a, 05.70.Fh, 64.60.Ak 



Quantum phase transitions (QPTs) have been a research 
topic subject to intense study, since their significant role was 
realized in accounting for high-T c superconductors, fractional 
quantum Hall liquids, and quantum magnets |HE]- Recently, 
significant advances have been made in attempt to clarify the 
connection between quantum many-body physics and quan- 
tum information science. This provides a newperspective 
to investigate QPTs from entanglement J2I H la, 0] and/z- 
delity B8l|9 Ul All ill 1211 . basic notions of quantum information 
science 11311 and turns out to be very insightful in our under- 
standing of QPTs in a variety of quantum lattice systems in 
condensed matter. 

Conventionally, orders and fluctuations provide a proper 
language to study QPTs, with order parameters being the key 
to quantify quantum fluctuations. Instead, the fidelity ap- 
proach is based on state distinguishability arising from the or- 
thogonality of different ground states in the thermodynamic 
limit. In fact, the ground state fidelity for a quantum sys- 
tem may be mapped onto the partition function of the equiv- 
alent classical statistical lattice model with the same geome- 
try ITlll . Thus, the fidelity per site is well-defined in the ther- 
modynamic limit, and its singularities unveil transition points, 
at which the system under consideration undergoes QPTs. 
Therefore, a practical means is now available to map out the 
ground state phase diagram for a quantum lattice system with- 
out prior knowledge of order parameters. An intriguing ques- 
tion is how to characterize singularities in the fidelity per site. 
Indeed, a proper answer to this question will shed new light 
on our understanding of QPTs. 

In this paper, we present an intrinsic characterization of sin- 
gularities in the fidelity per site in terms of Riemannian ge- 
ometry. For this purpose, we first define a fidelity surface as 
a surface embedded in a Euclidean space, which in turn is de- 
termined by the average fidelity per lattice site between two 
ground states of a quantum lattice system as a function of the 
control parameters. This makes the whole machinery devel- 
oped in differential geometry of curves and surfaces available 
to study QPTs. As it is well known, the Gaussian curvature, 
or equivalently, the Ricci scalar curvature for the surfaces em- 



bedded in Euclidean spaces, is a fundamental concept used to 
measure how curved a surface is. Therefore, the Gaussian cur- 
vature is expected to naturally quantifies quantum fluctuations 
that destroy orders at transition points. We discuss the global 
behaviors of the Gaussian curvature. It turns out that quantum 
fluctuations wildly distort the fidelity surfaces near the tran- 
sition points. Generically, precursors of QPTs occur in the 
Gaussian curvature for finite-size systems. In the thermody- 
namic limit, the Gaussian curvature becomes singular at tran- 
sition points. The one-dimensional quantum Ising model in a 
transverse field is exploited to explicitly illustrate the theory. 
We also perform a finite size scaling analysis for the Gaussian 
curvature with different lattice sizes to extract the correlation 
length critical exponent. 

Fidelity surfaces. For a quantum lattice system described 
by a Hamiltonian H(A), with A a control parameter. Here we 
restrict ourselves to discuss the simplest case with one single 
control parameter, although the extension to multiple control 
parameters is straightforward. For two ground states \tfr(Ai)) 
and \1ffiA2)) corresponding to different values of the control pa- 
rameter A, the fidelity is defined as F(Au^2) = \(iK^2)\^(Ai))l 
For a large but finite L, the fidelity F asymptotically scales as 
F(Ai,A2) ~ d L (A\,A2), where the scaling parameter d(A\,A 2 ) 
characterizes how fast the fidelity changes when the thermo- 
dynamic limit is approached ifioll . Physically, it is the fidelity 
per site. Here note that the contribution from each site to 
F(A\,A<£) is multiplicative. Following 111 111 , the ground state 
fidelity for a quantum system is nothing but the partition func- 
tion of the equivalent classical statistical lattice model with 
the same geometry, if one utilizes the tensor network repre- 
sentations of ground state many-body wave functions. There- 
fore, dfAi, A2) may be interpreted as the partition function per 
site 01411 . which is well-defined in the thermodynamic limit: 



\nd(A u A 2 ) = lim lnF(A u A 2 )/L. 



(1) 



The fidelity per site d{A\, A%) satisfies the properties: (1) sym- 
metry under interchange A \ < — > A 2 ; (2) d(A\,A\) = 1; and (3) 
< d(A u A 2 ) < 1. 

For simplicity, let us assume that the system undergoes a 



2 



QPT at A c . If \if/(Ai)} and \ip(A 2 )) are in the same phase, 
then they flow to the same stable fixed point in the sense of 
renormalization group theory, and so their difference arises 
from quantum fluctuations depending on the details of the 
system. On the other hand, if \ijj(A\)) and ^(Aq)) are in dif- 
ferent phases, then they flow to two different stable fixed 
points. Therefore, they possess different orders, although 
quantum fluctuations originate from the same unstable fixed 
point A c lfl5ll . Imagine that if there were no quantum fluc- 
tuations, then d{A\,A 2 ) would be simply 1 when \ifr(A\)) and 
\i//(A 2 )) are in the same phase; otherwise, when \iff(A\)) and 
\if/(A 2 )) are in different phases, d(A\,A2) would take the mini- 
mum value corresponding to the two stable fixed points. For 
continuous QPTs, quantum fluctuations are strong enough 
such that no orders survive at the transition point, so d( A\ , At) 
is continuous, but displays singularities, whereas for the first 
order QPTs, d{A\, A%) remains to be discontinuous at the tran- 
sition point. An interesting observation is to regard the fidelity 
per site, d{Ai,A 2 ), as a two-dimensional surface embedded 
in the three-dimensional Euclidean space, with a Riemannian 
metric induced from the Euclidean metric. Our aim is to give 
an intrinsic characterization of singularities in such a fidelity 
surface in terms of Riemannian geometry. 

Differential geometry of the two-dimensional surfaces em- 
bedded in the three-dimensional Euclidean space. Let us 
briefly recall the fundamentals of differential geometry of 
surfaces embedded in Euclidean spaces 111 611 . For a two- 
dimensional surface embedded in a three-dimensional Eu- 
clidean space: z = f(Ai,A 2 ), the first fundamental form on 
the surface is 



dl 2 = gij dA i dA j = E(duf + 2F(dudv) + G(dvf 



(2) 



where gy is the Riemannian metric on the surface: gn = 1 + 
Jf^gu = gn = fxjh, and g 22 = 1 + /£. Here the subscripts 
A\ and A% denote partial differentiations with respect to A\ and 
A2, respectively. In terms of the co-ordinates u - A\ and v = 
A 2 , we have E - g n , F = g 12 = gu and G = g 22 . Suppose 
the surface is given in parametric form: r — r(u, v). Then, the 
vector product r u x r v is a non-zero vector perpendicular to the 
surface at each non-singular point; define m to be a unit vector 
in the normal direction, then one has r u x r v = \r u x r v \m. For 
a curve r = r(u(l), v(l)) on the surface, the projection of the 
second order derivative r of r with respect to the arc length I 
on the normal to the surface leads to the second fundamental 
form as follows 

(r, m)(dlf = bijdA'dA j = X(duf + lYdudv + Z(dv) 2 , (3) 

if a surface is given in the form z = f(A\ , A 2 ) with A\ — u, A 2 - 
v, and r(u,v) - (u, v,f(u, v)). Therefore, we have X — b\\ - 

/WV*+4 + fl> Y = bn = b 2l = f AtAl /^l + fl+fl 

and Z = b 22 = f hA J ^jl+fl +ft_. 

The eigenvalues of the pair of quadratic forms (f2]l and (01 
are the principal curvatures of the surface at the point under 
investigation. The product of the principal curvatures is the 



Gaussian curvature K of the surface at the point, and their 
sum the mean curvature. The principal curvatures k\ and k 2 
are the solutions of equation: 



det(2 - kG) = 0, 



(4) 



where Q = (foy) is the matrix of the second fundamental form, 
and G = igij). Since the first fundamental form is positive 
definite, its matrix G is non-singular. Hence det(Q - kG) = 
detG det(G l Q -k- 1), we deduce that the Gaussian curvature 
K — k\k 2 — det(G~'0 = detQ/detG and the mean curvature 
M = k\ + k 2 — tx(G l Q). Therefore, the Gaussian curvature K 
and the mean curvature M take the form: 



K = 



fX\X\fXi\i f 



■2 

A.1A.2 



and 



M 



( 1+ fl + fl) 



(5) 



(6) 



respectively. We notice that the sign of the Gaussian curva- 
ture K is the same as the sign of the determinant: f^Xif^M ~ 
ff lAa , i.e., the Hessian of z — f(A\ , A 2 ) 

It follows that, in contrast with the mean curvature M, the 
Gaussian curvature K of a surface may be expressed in terms 
of the induced metric on the surface alone, and is therefore 
an intrinsic invariant of the surface In addition, a two- 
dimensional surface in a three-dimensional space may also be 
regarded as a differentiable manifold endowed with a Rieman- 
nian metric induced from the Euclidean metric. The Ricci 
scalar curvature R is twice the Gaussian curvature K: R — 2K. 

Global behaviors of the Gaussian curvature K for a fi- 
delity surface. Now we consider the (logarithmic function 
of) fidelity per site, In d(A\ , A 2 ), as a two-dimensional sur- 
face embedded in the three-dimensional Euclidean space: z = 
f(A\,A 2 ) = lnd(A\,A 2 ). The Gaussian curvature K(A\, A 2 ) 
for such a fidelity surface may be used to quantify how 
strong quantum fluctuations are in given quantum many-body 
ground states, thus providing an intrinsic characterization of 
singularities in the fidelity surface. Indeed, as justified in 
Refs. $M 

1 111 , the fidelity per site d(A\,A 2 ) is singular when 
A\{A 2 ) crosses A c for a fixed A 2 (Ai) in the thermodynamic 
limit. Therefore the Gaussian curvature K(A\,A 2 ) for the fi- 
delity surface is singular at A\ — A c and/or A 2 = A c in the ther- 
modynamic limit. Generically, we have: (1) K(A\,A 2 ) > 0, 
there is a neighborhood of the point throughout which the sur- 
face lies on one sides of the tangent plane at the points; (2) 
K(Ai,A 2 ) < 0, then the surface intersects the tangent plane at 
the point arbitrarily close to the point. If the surface is strictly 
convex, then we say that the Gaussian curvature K{A\,A 2 ) is 
positive at every point of the surface. That is what happens 
if A\ and A 2 are away from the transition point. However, if 
A\ and A 2 are close to the transition point, then the Gaussian 
curvature K(A\,A 2 ) can be negative. 
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For finite-size systems, the Gaussian curvature K(Ai , A2) re- 
mains to be smooth, although the precursors of QPTs occur as 
anomalies in the Gaussian curvature K(A 1 ,A2). The anoma- 
lies get more pronounced when the thermodynamic limit is 
approached. We may take advantage of this fact to perform 
finite size scaling to extract the correlation length critical ex- 
ponent. 




FIG. 1 : (color online) The behavior near the critical point A c = 1 is 
analyzed for the Gaussian curvature K(A 1 ,A 2 ) of the quantum trans- 
verse Ising model for various lattice sizes. The curves shown cor- 
respond to different lattice sizes L = 201,401, 1201,2001, and 00. 
The peaks (dips) get more pronounced in the left (right) side with 
increasing system size. The Gaussian curvature K(A t , A 2 ) diverges at 
the critical point A t = A c for the infinite-size system (L = 00). Upper 
panel: Here K(A t , A 2 ) is regarded as a function of A t for A 2 = 0.6 and 
y=l. Lower panel: Here K(Ai,A 2 ) is regarded as a function of A\ 
for A 2 = 0.6 and y = 1/2. 

Quantum XY spin 1/2 model. The quantum XY spin model 
is described by the Hamiltonian 



M 



1 + y 



1 



(7) 



Here <j x , cr. and cr" are the Pauli matrices at the ;'-th lattice 

j' j j J 

site. The parameter y denotes an anisotropy in the nearest- 
neighbor spin-spin interaction, whereas A is an external mag- 
netic field. The Hamiltonian (0 may be exactly diagonal- 
ized 11171 ll 811 for any finite size L with L = 2M + 1 . In the 
thermodynamic limit L — > 00, In d(A\ , A2) takes the form Q : 



1 r" 

lnd(A[,A2) = — I da In T (A\, A2; a), 
2x Jo 



(8) 



where T(A\,A2\d) = cos[)?(/l 1 ; a) - ^{A^, a)]/2, with 
cos i?(/l; a) = (cos a - A)/ -^(cos a - A) 2 + y 2 sin 2 a 1 191. 



Now it is straightforward to calculate the Gaussian cur- 
vature K(A\,A2) for the fidelity surface of the quantum XY 
spin chain. In Fig. [T] we plot K(A\,A2 = 0.6) for the fi- 
delity surface of the quantum XY model (y = 1 for the upper 
panel and y — 1/2 for the lower panel). One observes that 
K(Ai,A2 = 0.6) is divergent as a function of A\ at the criti- 
cal point A c = 1 for the infinite-size system L — 00, indicat- 
ing that the fidelity surface is wildly distorted, due to strong 
quantum fluctuations near the critical point. This is true for 
any nonzero y, consistent with the fact that the quantum XY 
model for any nonzero y belongs to the same universality class 
as the quantum transverse Ising model. That is, there is a crit- 
ical line y + and A c — 1 ; only one (second-order) critical 
point Ac - 1 separates two gapful phases: (spin reversal) Z2 
symmetry-breaking and symmetric phases. 
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FIG. 2: (color online) (a) The peaks values of the Gaussian curvature 
K(Ai , A 2 ) of the quantum transverse Ising model for large lattice sizes 
scale as L/(lnL) 4 . (b) The dips values of the Gaussian curvature 
K(A\ , A 2 ) of the quantum transverse Ising model for large lattice sizes 
scale as L/(lnL) 4 . In both cases, A 2 = 0.6 and y = 1. 

Finite size scaling analysis for the Gaussian curvature K. 
We focus on the quantum Ising universality class. The order 
parameter, i.e., magnetization (cr x ) is non-zero for A < 1, and 
otherwise zero. At the critical point, the correlation length 
f ~ \A — /t c r v with v = 1 11811 , In order to analyze how the 
Gaussian curvature ^(^1,^2) behaves near the critical point 
A c — 1, we perform a finite size scaling analysis for the quan- 
tum transverse Ising model. 

As already observed, the drastic change of the ground state 
wave functions makes the Gaussian curvature K(A\, A2) diver- 
gent when the system undergoes the second order QPT at the 
critical point A c — 1 in the thermodynamic limit. However, 
for finite-size systems, K(Ai,A2) remains to be smooth for 
the quantum XY model. In Fig. [TJ the numerical results are 
also plotted for the Gaussian curvature K(A\,A2) with differ- 
ent system sizes, where A2 = 0.6 and y = 1 (upper panel) and 
A2 = 0.6 and y = 1/2 (lower panel). More precisely, in the 
thermodynamic limit, K(Ai , A2) (as a function of Ai for a fixed 
A2) diverges at the critical point A\ — A c : 



K(A U A 2 ) 



1 



(9) 



l/l! -/LJ(ln 1^-^1)4 • 

However, there is no divergence for finite-size systems, but 
there are clear anomalies, featuring two quasi-critical values 
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A p and Ad, one at each side of the critical point. On the 
left (right) side, the so-called quasi-critical points A p (Ad) ap- 
proach the critical value as A p « 1 - 1.6149L _103531 (/l d « 
1 + 9.69198ZT 974152 ), with the values at peaks (dips) diverg- 
ing with increasing system size L: 



K(A u A 2 )\=k 



(InL) 4 



+ constant. 



(10) 



Here the prefactor k p (d) is non-universal in the sense that it 
depends on A 2 and y. We emphasize that Eq. ( fTOb follows 
from Eq. {9), if we take into account the fact that the model 
is conformally invariant at the critical point. Indeed, on the 
one hand, from Eq. © and the correlation length £ ~ \A. — 
A c \- V with v = 1, we have K{A X ,A 2 ) ~ <f/(ln£) 4 . On the other 
hand, the conformal invariance requires the scale invariance: 
g/L = i; jU . The numerical results are, respectively, plotted 
*-m m Fig. |2]and for A p ^) in Fig. [3] with A 2 - 
0.6 and y — 1. The same is also true for any nonzero y. This 
shows that, consistent with the exact result, the correlation 
length critical exponent v equals 1, as long as y is nonzero. 





FIG. 3: (color online) (a): The positions of the peaks approach 
the critical point A c = 1 with increasing system size L as A p at 
1 - 1.61490L" 1 03531 . (b) The positions of the dips approach the 
critical point A c = 1 with increasing system size L as Aj at 1 + 
9.69198ZT 974152 . In both cases, A 2 = 0.6 and y = 1. 

Conclusions. We have shown that singularities in fidelity 
surfaces may be intrinsically characterized in terms of Rie- 
mannian geometry, based on the fidelity description of QPTs. 
Genetically, the Ricci curvature tensor for finite-size systems 
is analytic and it exhibits singularities at transition points 
in the thermodynamic limit, as reflected in the Ricci scalar 
curvature that blows up when the system size tends to oo. 
This opens up the possibility to exploit the theory of Ricci 
flows 12011 to characterize QPTs in condensed matter theory. 
The one-dimensional quantum Ising model in a transverse 
field is exploited as an example to explicitly illustrate the the- 
ory 12111 . and a finite size scaling analysis has been performed 
for the Ricci scalar curvature with different lattice sizes, and 
the correlation length critical exponent has been extracted, 
consistent with the known exact value. 
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